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ABSTRACT 

A  fully  coupled  mathematical  model  describing  the  interactions  between  a  vibrating  thin 
cylindrical  shell  and  an  enclosed  acoustic  field  is  presented.  Because  the  model  will  ultimately 
be  used  in  control  applications  involving  piezoceramic  actuators,  the  loads  and  material  contri¬ 
butions  resulting  from  piezoceramic  patches  bonded  to  the  shell  are  included  in  the  discussion. 
Theoretical  and  computational  issues  lead  to  the  consideration  of  a  weak  form  of  the  model¬ 
ing  set  of  partial  differential  equations  (PDE’s)  and  through  the  use  a  semigroup  formulation, 
well-posedness  results  for  the  system  model  are  obtained. 


’The  research  of  H.T.B.  was  supported  in  part  by  the  Air  Force  Office  of  .Scientific  Research  under  grant 
AFOSR-90-0091.  This  research  was  also  supported  by  the  National  Aeronautics  and  Space  Administration 
under  NASA  Contract  Numbers  NAS1-I860.5  and  NASI- 19480  while  H.T.B.  was  a  visiting  scientist  and 
R.C.S.  was  in  residence  at  the  Institute  for  Computer  Applications  in  Science  and  Engineering  (ICASE), 
NASA  Langley  Research  Center,  Hampton,  VA  2.3681. 


1  Introduction 


The  active  control  of  noise  generated  by  structural  vibrations  has  been  studied  for  several 
years  in  various  applications.  One  recent  motivation  leading  to  an  intense  study  of  prob¬ 
lems  involving  the  reduction  of  structure-borne  noise  has  resulted  from  the  development 
of  a  new  class  of  turboprop  and  turbofan  engines.  These  engines,  although  very  fuel  ef¬ 
ficient,  are  also  very  noisy.  Specifically,  the  low  frequency,  high  magnitude  exterior  noise 
fields  produced  by  these  new  engines  cause  vibrations  in  the  fuselage  which,  due  to  struc¬ 
tural  acoustic  coupling  between  the  fuselage  and  the  interior  acoustic  field,  lead  to  un¬ 
acceptably  high  cabin  noise  levels.  The  problem  is  exacerbated  by  the  increased  use  of 
lightweight,  composite  materials  in  the  cabin  walls.  In  order  to  make  the  use  of  the  en¬ 
gines  feasible  in  commercial  aircraft,  a  great  deal  of  research  has  been  aimed  at  developing 
control  techniques  to  reduce  this  unwanted  interior  noise. 

Control  techniques  for  structural  acoustics  problems  of  this  type  have  been  studied  from  a 
variety  of  perspectives  [2,  7,  12,  13,  14,  15,  16,  20,  21,  25],  one  of  which  is  through  the  use  of 
piezoceramic  patches  which  are  bonded  to  the  enclosing  structure  [2,  7,  14,  16).  As  detailed 
in  [8,  9],  the  patches  create  bending  moments  and  in-plane  strains  when  a  voltage  is  applied, 
and  through  these  actions,  can  be  used  to  alter  the  structural  dynamics  in  a  manner  which 
ultimately  reduces  interior  noise. 

In  this  paper,  we  develop  a  fully  coupled  mathematical  model,  as  well  as  well-posedness 
results  for  this  model,  which  can  be  used  when  applying  parameter  estimation  and  PDE^based 
control  strategies  to  problems  involving  the  use  of  piezocc  ^mic  actuators.  This  model  extends 
the  2-D  results  in  [2,  7]  and  the  3-D  results  in  [10]  to  the  cylindrical  domain  of  interest  (the 
3-D  domain  in  [10]  consisted  of  a  hard-walled  cavity  with  a  thin  plate  at  one  end).  While 
ultimately  motivated  by  the  cabin  noise  problem  mentioned  above,  the  model  is  designed  to 
be  consistent  with  an  experimental  apparatus  being  used  at  the  Acoustic  Division,  NASA 
Langley  Research  Center.  This  setup  consists  of  a  thin-walled  hollow  cylindrical  shell  which 
is  supported  by  rigid  caps  at  the  ends.  Piezoceramic  patches  bonded  to  the  walls  of  the  shell 
are  then  used  to  control  interior  noise  which  has  been  generated  by  the  vibrations  of  the  walls. 

This  mathematical  model  differs  from  previously  used  3-D  shell/acoustic  models  in  that 
it  fully  incorporates  the  backpressure  and  momentum  conditions  which  couple  the  structural 
dynamics  with  the  interior  acoustic  response,  thus  yielding  a  time-dependent  system  of  partial 
differential  equations  which  describe  the  system  dynamics.  Because  the  acoustic/structure 
interactions  are  fully  incorporated,  this  model  is  useful  for  control  strategies  which  utilize  the 
natural  “feedback”  loop  due  to  the  coupling  between  the  acoustic  fields  and  the  structural 
vibrations.  Moreover,  by  incorporating  the  piezoceramic  patch/shell  interaction  results  of  [8], 
the  model  can  be  employed  in  devising  control  techniques  utilizing  piezoceramic  actuators. 

After  the  development  of  the  mathematical  model,  well-posedness  results  are  presented. 
These  are  important  not  only  for  determining  the  existence  and  continuous  dependence  of 
solutions,  but  also  in  providing  an  initial  framework  that  can  be  used  when  determining 
suitable  approximation  schemes.  By  carefully  noting  the  Hilbert  spaces  containing  the  state 
variables  and  test  functions,  appropriate  choices  can  be  made  for  the  approximating  subspares 
which  contain  the  basis  functions  used  in  finite  element  or  spectral  approximations.  This  is 
important  not  only  when  considering  forward  simulations  but  also  when  considering  theoretical 
issues  concerning  the  parameter  estimation  and  control  problems. 


Section  2  of  this  presentation  contains  a  description  of  the  components  of  this  system 
along  with  the  strong  form  of  the  equations  of  motion  for  the  system.  The  disadvantages  of 
the  strong  form  are  discussed  in  Section  3  and  a  weak  or  variational  form  of  the  system  equa¬ 
tions  is  presented.  An  abstract  framework  amenable  to  forward  simulations,  the  estimation 
of  physical  parameters,  and  the  implementation  of  PDE-based  feedback  control  strategies  is 
also  developed  in  this  section.  In  Section  4,  the  abstract  model  described  in  the  third  section 
is  shown  to  be  well-posed.  This  is  accomplished  using  “extrapolation  space”  ideas  and  argu¬ 
ments  similar  to  those  presented  in  [5,  6,  17].  Having  obtained  the  existence  and  continuous 
dependence  of  solutions  for  the  model,  approximation  and  LQR  optimal  control  techniques 
similar  to  those  discussed  in  [2,  7,  10]  can  be  applied  to  the  problem  in  an  effort  to  suitably 
reduce  interior  pressure  levels. 


2  Strong  Form  of  the  System  Equations 

As  motivated  by  the  experimental  setup  described  in  the  Introduction,  the  structural  acoustics 
problem  under  consideration  is  assumed  to  consist  of  a  cylindrical  acoustic  cavity  Q{t)  which 
is  enclosed  by  a  thin  cylindrical  shell  (see  Figure  1).  Hard  wall  conditions  are  taken  at  the 
ends  r  of  the  cavity  in  order  to  model  the  rigid  end  caps  used  in  the  experimental  apparatus. 
These  end  caps  also  dictate  the  use  of  clamped  boundary  conditions  when  approximating  the 
shell  dynamics. 

As  outlined  in  the  discussion  in  the  previous  section,  the  dynamics  of  the  coupled  system 
are  composed  of  an  acoustic  response,  shell  dynamics,  piezoceramic  patch/shell  interactions, 
and  the  coupling  acoustic/structure  interactions.  Each  of  these  components  is  briefly  described 
and  the  results  are  summarized  in  a  coupled  system  of  PDE’s  which  describe  the  system 
dynamics. 


Figure  1.  Cylindrical  acoustic  cavity  H{t)  with  hard  end  caps  F. 

2.1.  Acoustic  response 

The  acoustic  wave  motion  inside  the  cavity  can  be  described  either  in  terms  of  a  velocity 
potential  (j)  or  the  acoustic  pressure  p  (the  two  are  related  through  the  relationship  p  = 
P](i>i  where  pj  is  the  equilibrium  density  of  the  atmosphere);  motivated  by  control  theoretic 
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considerations  and  to  simplify  the  presentation  which  follows,  we  choose  the  former  as  the 
second-order  state  variable.  The  acoustic  dynamics  inside  the  cavity  are  then  modeled  by  t’.e 
undamped  wave  equation 

<l>tt  =  ,  (r,0,x)  e  >  0  , 

V(f>'h  =  0  ,  {r,6,x)  E  r  ,t  >  0 

where  c  is  the  speed  of  sound  in  the  cavity,  h  is  the  outward  axial  unit  normal  to  the  cavity 
end  caps,  and  the  Laplacian  in  cylindrical  coordinates  is  given  by 

.  .  d^(l>  1  d(l>  1  d'^(l)  d'^(f) 

^  dr^  r  dr  80^  dx^ 

Again,  P  denotes  the  circular  caps  at  the  ends  of  the  cavity  Q{t)  (the  cavity  is  variable  in 
time  due  to  the  fact  that  the  enclosing  shell  is  vibrating),  and  x  is  the  axial  direction  in  the 
cavity.  Finally,  we  note  that  air  damping  inside  the  cavity  was  omitted  due  to  the  relatively 
small  dimensions  of  the  experimental  cavity,  and  could  readily  be  incorporated  for  significantly 
larger  cavities. 

2.2.  Shell  dynamics 

In  order  to  specify  the  equations  of  motion  of  the  bounding  structure,  we  consider  a  thin 
cylindrical  shell  having  length  thickness  h  and  radius  R  with  axial  direction  taken  along 
the  x-axis  (see  Figure  2).  The  displacements  of  the  middle  surface  in  the  axial,  tangential  and 
radial  directions  are  taken  to  be  u,  v  and  w,  respectively.  Furthermore,  the  shell  is  assumed  to 
have  density  p.  Young’s  modulus  E,  Poisson  ration  u,  and  damping  coefficient  Cp-  Finally,  we 
assume  that  the  shell’s  length  is  relatively  short  in  relation  to  it  radius  and  hence  the  Donnell- 
Mushtari  equations  can  be  used  when  approximating  its  motion  [18].  This  last  assumption 
is  made  purely  for  ease  of  presentation  and  higher-order  theories  such  as  the  Byrne- Fliigge- 
Lur’ye  model  can  be  substituted  for  that  of  Donnell  and  Mushtari  as  warranted  by  the  shell 
dimensions. 


□  □ 


As  presented  in  [8,  19],  the  Donnell-Mushtari  equations  for  a  thin  cylindrical  shell  are  given 


„  d'^v  dNe 

-  Se 


D  ^  ^  ^  ^  O 


dNg, 

de 

dx 

d{Ng)pe  ^ 

^  dx 

de  ‘  ’ 

d^Mxg 

+  Ng  =  Rqn  — 

dxdO 

S,,2{x)‘k2{0) 


s,, 2(^)3, M 


,d'^{M:,)pe 
'  dx^ 


1  d\Me)p. 
R  de^ 


Here  A/i,  A/g,  and  Ng  are  internal  moments  and  force  resultants,  (A/i)pe,  {Mg)pe,  (A^x)pe 
and  {Ng)pe  are  the  resultants  for  the  loads  generated  by  the  piezoceramic  patches  when  a 
voltage  is  applied,  and  %  is  the  external,  normal  load  on  the  shell.  For  a  patch  with  uniform 
thickness  and  bounding  values  Xi,X2,0^  and  $2  as  depicted  in  Figure  2,  the  presence  of  the 
indicator  function 

'1  ,  X  <  (li  +  X2)/2 

S\,2ix)  =  -  0  ,  X  =  (xi  +  X-2)/2  (2) 

-1  ,  X  >  (xi  +  X2)/2 

derives  from  the  fact  that  the  forces  generated  by  the  patch  in  the  x-direction  are  antisym¬ 
metric  (equal  in  magnitude  but  opposite  in  sign)  about  the  line  x  =  (xi  -|-  X2)/2.  The  same 
holds  true  for  the  forces  in  the  0-direction  and  6V2(0)  is  defined  in  an  analogous  manner. 

From  the  discussion  in  [9],  internal  damping  can  be  incorporated  in  the  shell  equations  by 
assuming  a  constitutive  law  which  posits  that  stress  is  proportional  to  a  linear  combination  of 
strain  and  strain  rate.  This  yields  a  a  Kelvin-Voigt  type  of  damping  in  the  shell  with  internal 
resultants  for  shell  regions  not  covered  by  patches  given  by 


Eh  [cJtx  f  \  dv  w\ 
"  (1  -i/^)  ■*'  RJ 


Cph  d  du  ( \  dv  w 

(1  -1/2)^  [r'^  ^  R 


Eh  r  1  w  du\  Cnh  f  1  w  du 
®  (1-1/2)  dx\  {\-u^)dt[Rde  R  dx 

Eh  dv  \  du  ^  C[)h  d  dv  \  du 

N^e  =  Ng^  =  +  2(1  +  ,2)^ 

Eh^  d^w  u  d^w  Cph^  d  d^w  u  d^w 

^  “12(1  -  //2)  ~  12(1  -7/2)^  WW 

Eh^  1  d^w  d^w}  Cph^  d  1  d^w  d^tv 
""  “l2(l  -7/2)  [WW  “  12(1  -1/2)^ 

Eh^  d^w  Cph^  d  d^w 
M^g  =Mex  =  ~  i2Ri^\  ^  y)  dxdd  ~  V2R.{l+u)Ft  [dlde  ■ 
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Similar  expressions  (e.g.,  see  [9])  can  be  derived  for  those  regions  of  the  structure  in  which 
piezoceramic  patches  are  bonded  to  the  shell.  In  those  areas,  the  internal  force  and  mo¬ 
ment  resultants  contain  contributions  not  only  from  the  shell  but  also  from  the  internal 
stresses  in  the  piezoceramic  patches.  The  contributions  resulting  from  the  external  loads 
{Mx)pe,  {M8)pe,  {Nx)pe,  {No)pe  and  q,i  will  be  discussed  further  in  the  next  two  subsections. 

Finally,  we  note  that  the  boundary  conditions  for  the  shell  must  be  consistent  with  the 
rigid  end  cap  conditions  on  the  experimental  apparatus.  Depending  on  the  exact  nature  of 
the  end  caps,  one  may  enforce  either  the  clamped  edge  condition 

dtv 

u  =  V  —  w  —  =  U 

ox 

or  the  simply  supported  edge  condition 

u  =  V  =  w  =  Mx  =  0 

at  the  ends  of  the  shell.  For  definiteness,  we  will  employ  the  first  while  simply  noting  that 
the  latter  often  adequately  models  the  conditions  being  maintained  at  the  edges  of  the  shell. 


2.3.  Piezoceramic  patch/cylindrical  shell  interactions 

In  order  to  specify  the  moment  and  force  resultants  which  result  from  the  activation  of  the 
piezoceramic  patches,  we  consider  patches  which  are  bonded  to  the  inner  and  outer  surfaces 
of  the  shell  in  a  manner  such  that  their  edges  are  parallel  to  lines  of  constant  x  and  0  as  shown 
in  Figure  2.  When  it  is  necessary  to  differentiate  between  the  two  patches,  the  outer  will  be 
denoted  with  a  subscript  pei  with  a  subscript  pe-z  being  used  to  denote  the  inner  patch  in 
each  pair.  We  point  out  that  in  order  to  keep  this  presentation  compact,  we  have  assumed 
that  both  patches  are  active,  and  we  refer  the  reader  to  [8]  for  the  case  when  only  one  patch 
is  bonded  to  the  shell. 

As  discussed  in  [8],  the  total  line  moments  and  forces  for  a  single  patch  pair  with  bounding 
values  and  02  are 

{Mx)p,  =  [{MoU  +  iMe)p.A  [H,{x)  -  Hzix)]  [H,{0)  -  Hz{0)] 
iMe)pe  =  [{Mx)pe,  +  iMx)pe,]  [H^{x)  -  Hz{x)]  [H,{0)  -  Hz{0)] 

(4) 

{Nx)pe  =  [{Nx)pe,  +  {Nx)pe,]  [H,{x)  -  Hzix)]  [H,{0)  -  Hz(0)]  S,,z{x)S^,z{0) 


{Ne)pe  =  [(A^0)pe.  +  (A^^)peJ  {^,{0)  -  S,Ax)Su2{0)  . 

Here  H  is  the  Heaviside  function  with  Hi{x)  =  H{x  —  a:,),i  =  1,2,  and  .S’,  2  is  the  indicator 
function  defined  in  (2)  with  similar  definitions  being  used  in  0.  The  individual  components 


are  given  by 


i^e)pe\  —  1  _  ,.2  *''pe)Cpe, 

'■  ‘^pe 


{^e)pe2  —  1  _  ,.2  *v)^p«’l  • 

where  Epe,  t'pe  and  T  are  the  Young’s  modulus,  Poisson  ratio  and  thickness  of  the  patches, 
respectively  (it  is  assumed  that  the  inner  and  outer  patches  have  the  same  material  properties). 
With  ^31,  Vi  and  V-i  denoting  the  piezoceramic  strain  constant  and  voltages  into  the  two 
patches,  the  terms  Cpe,  =  ^V\  and  Cpej  =  ^V-z  provide  a  relationship  between  the  applied 
voltages  and  the  resulting  in-plane  mechanical  strains. 

From  (4)  and  (5),  it  can  be  seen  that  the  resultants  due  to  the  activation  of  the  patches 
depend  on  the  material  and  geometric  properties  of  the  patches,  the  radius  of  the  shell,  and 
the  voltage  being  applied  to  the  patches.  Once  determined,  these  resultants  are  substituted 
into  (1)  as  the  external  loads  on  the  shell.  If  multiple  patch  pairs  are  present,  the  resultants 
for  each  pair  are  incorporated  in  the  shell  equations  in  a  similar  manner. 

2.4.  Normal  loads  and  coupling  conditions 

It  can  be  seen  in  (1)  that  the  shell  equations  also  contain  the  term  which  represents  the 
normal  load  on  the  shell.  In  writing  the  equations  in  this  form,  we  have  made  the  assumption 
that  the  only  loads  on  the  shell  are  the  contributions  due  to  the  patches,  and  the  normal 
load  which  is  due  to  the  exterior  noise  field  /  as  well  as  backpressure  generated  by  the  interior 
acoustic  field.  This  assumption  was  made  merely  to  simplify  the  discussion,  and  more  comi)lex 


acoustic  or  mechanical  loads  on  the  shell  can  be  incorporated  as  additional  external  forces  and 
moments  in  the  shell  equations.  Because  in  this  case,  consists  of  contributions  from  the 
exterior  and  interior  acoustic  fields,  the  first  coupling  relation  is 

•^)  —  y(^)  Pf^l  (^1  3r),  v{t^  9^  jj),  u{t^  0,  j)) 

and  is  in  general  nonlinear  since  the  backpressure  onto  the  shell  occurs  at  the  shell’s  surface. 
The  second  coupling  condition  is  the  velocity  constraint 

^  (t,u;(t,0,a;),u(t,0,  x),u(t,0,a!:))  =  Wt{t,9,  x) 
or 

which  simply  states  that  the  shell  is  impermeable  to  air.  Note  that  this  second  constraint  also 
provides  a  boundary  condition  for  the  acoustic  response. 

2.5.  Fully  coupled  system:  strong  form 

Under  the  assumption  of  small  displacements  which  is  inherent  in  the  linear  shell  theories, 
the  variable  domain  fl(f)  is  approximated  by  the  fixed  domain  fl,  and  the  general  nonlinear 
coupling  conditions  are  approximated  by  their  linear  approximations.  Let  the  boundary  To 
denote  the  shell  in  its  unperturbed  state.  For  s  pairs  of  piezoceramic  patches,  this  then  yields 
the  linear  model 

=  ,  {r,9,x)  e  0.  ,t  >  0  , 

\7<f>-n  =  0  ,  (r, 0, x)  G  r  >  0 


^{t,R,e,x)  =  Wt{t,6,x)  ,  (0,x)€ro,t>O 


d'^U  dNer  ny^d[iN^)pe]i 


^  dx 


69 


=  -«E 


i=l 


dx 


1=1 

„  ,  s^w  1  d^M,  „d‘Ma  ,  „ 

=  R  [f{t,  9,  x)  -  p/<A{(t,  R,  9,  x)] 


1=1  I 


/  1  d^jMeU, 

dx^  R?  69^ 


dw 

u  =  V  =  w  =  — —  =  0 
dx 


X  =  0,  f 


(l>(0,r,9,x)  =  (f>o(r,9,x)  ,  <?i((0,r,^,x)  =  <|i,(r,<9,x) 

u{0,9,x)  =  uo{9,x)  ,  u<(O,0.x)  =  ui(^,x) 

tj(O,0,x)  =  vo(^,a;)  ,  i;t(O,0.x)  =  t;i(0,x) 

«;(O,0,x)  =  tco(^,  a:)  ,  Wi{0,9,x)  =  wi{9,x) 


,  (r,9,x)efl 

{9,x)ero  . 


(0,x)  e  To  ,t  >  0 


(6) 
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Recall  that  the  internal  shell  moments  and  forces  are  summarized  in  (3)  while  the  expressions 
for  the  external  loads  generated  by  the  pair  of  piezoceramic  patches  are  given  in  (4).  The 
notation  [5'],2(^)]t  and  denotes  the  indicator  functions  (see  (2))  which  are  centered 

over  the  pair  of  patches.  Note  that  this  representation  admits  different  voltages  and 
geometries  for  the  .s  patches,  thus  increasing  the  flexibility  of  the  model  for  control  purposes. 

We  also  reiterate  that  the  only  damping  in  the  system  is  the  internal  Kelvin- Voigt  damping 
in  the  shell  although  medium  damping  inside  the  cavity  can  be  added  if  cavity  dimensions 
become  significant.  Without  medium  damping,  however,  the  system  is  only  weakly  damped, 
and  this  must  be  considered  when  discussing  theoretical  convergence  and  well-posedness  results 
for  the  problem. 


3  Weak  Form  of  the  System  Equations 

As  can  be  seen  in  the  model  (6),  the  use  of  the  strong  form  of  the  system  equations  leads  to 
first  and  second  derivatives  of  both  the  internal  and  external  moment  and  force  resultants. 
This  leads  to  difficulties  both  in  approximating  the  behavior  of  the  system  and  in  solving  the 
control  problem.  The  first  problem  results  from  the  presence  and  differing  material  properties 
of  the  piezoceramic  patches.  As  noted  in  [11]  where  piezoceramic  patches  are  bonded  to 
a  beam,  the  material  parameters  of  the  combined  structure  must  be  modeled  as  piecewise 
constants  in  order  to  accurately  match  structural  frequencies.  Hence  the  parameters  p,  £, 
and  Co  are  expressed  in  terms  of  a  Heaviside  basis  with  the  edges  of  the  patches  defining  the 
support  of  the  functions.  This  leads  to  difficulties  in  the  strong  form  since  it  necessitates  the 
differentiation  of  discontinuous  material  parameters.  A  similar  proble!Ti  arises  when  including 
the  moments  generated  by  the  patches.  As  seen  in  (4),  the  support  of  the  contributions  is 
given  in  terms  of  Heaviside  functions  which  implies  that  the  use  of  the  strong  form  yields  an 
unbounded  control  operator  (it  involves  the  differentiation  of  the  Heaviside  functions  as  well 
as  the  Dirac  delta).  To  avoid  these  difficulties,  it  is  advantageous  to  formulate  the  problem  in 
weak  or  variational  form  (the  use  of  the  variational  form  also  admits  the  use  of  basis  functions 
having  less  smoothness  than  those  used  when  approximating  the  solution  to  the  strong  form 
of  the  equations). 

3.1.  Weak  formulation 

The  second-order  state  for  the  problem  is  taken  to  be  y  =  (0,  u,  u,  u;)  in  the  Hilbert  space 
H  =  Z^(D)  X  L^(ro)  X  L^(ro)  X  L^(ro).  The  choice  of  the  space  Z^(D),  defined  as  the  quotient 
of  over  the  constant  functions,  results  from  the  fact  that  the  potentials  are  determined 

only  up  to  a  constant. 

To  provide  a  class  of  functions  which  are  considered  when  defining  a  variational  form  of 
the  problem,  we  also  define  the  Hilbert  space  V  =  x  //^(Tq)  x  //^(Tq)  x  //^(Fq)  where 

is  the  quotient  space  of  over  the  constant  functions.  The  subscript  0  in  the 

remaining  components  of  the  product  space  denotes  the  subset  of  functions  in  the  traditional 
Sobolev  spaces  which  satisfy  the  essential  boundary  conditions  u  =  u  =  u;  =  ^  =  0  at 
X  =  0,L 
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A  complete  discussion  concerning  the  derivation  of  the  a  variational  formulation  of  the 
shell  equations  from  energy  principles  can  be  found  in  [8].  For  our  purposes  here,  we  simply 
note  that  integration  in  combination  with  the  use  of  Green’s  theorem  yields  the  second-order 
variational  system 


f  I  ,  _  1  _  d'^m  1  2  ,  _  1  , 

■"1 1 -  f"* W  - 

^  '  (7) 

=  t  {[(^.)..i,  t  i  KW.  ^  -  i(«.).i,  §  -  ^  mu  ^ 

+  /  fvsdy 

Jro 

for  all  {^•,rii,T]2,T}z)  €  V  (note  that  du)  =  r dr dO dx  and  ^7  =  Rd9dx).  The  complex  inner 
products  are  used  in  anticipation  of  the  possibility  of  employing  complex  Fourier  expansions 
in  0  when  developing  suitable  approximation  schemes  for  the  problem. 

We  note  that,  as  is  usual,  in  this  variational  form,  the  derivatives  have  been  transferred 
from  the  plate  and  patch  moments  onto  the  test  functions.  This  eliminates  the  problem  of 
having  to  approximate  the  derivatives  of  the  Heaviside  function  and  the  Dirac  delta  which 
arises  in  the  case  of  the  strong  form  of  the  equations. 

The  system  (7)  can  be  formally  approximated  by  replacing  the  state  variables  by  their 
finite  dimensional  approximations  and  constructing  the  resulting  matrix  system.  In  order  to 
discuss  convergence  results  for  the  approximation,  parameter  estimation,  and  control  prob¬ 
lems,  however,  it  is  advantageous  to  pose  the  problem  in  terms  of  sesquilinear  forms  and  the 
bounded  operators  which  they  define. 


dx^  ”39'^  R  ^”3x39 


3.2.  Abstract  weak  formulation 

Before  proceeding  with  the  abstract  formulation,  it  is  necessary  to  describe  the  inner 
products  for  the  spaces  H  and  V  since  these  will  be  needed  when  determining  the  continuity 
of  various  operators.  From  energy  considerations,  it  is  appropriate  to  use  the  inner  products 

—  I  +  I  phurfid'y  +  f  phvlj^d'y  +  I  phwrj^d'f 

c‘‘  2  To  2  To  Jr  a 
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and 


=  f  Pf'Vcp-'V^dijJ 
Jn 

f  Eh  {\du  u  fdv  \1  di]\  1  ,,  ,■ 

ir.  (T^  life  "vj  2fl‘  . 


dv  1  du1  di]\  1 
dx  ^  R  06  _  06  j 


^  f  Eh  \  \  \  \  Ov  w  ^  ..^“1  dj)2  ,  1^,  .  X  1  diji]  , 

Ju  (1  1  [r'^  ^  ^  [^^Rm 


f  Eh  J  1  r  1  du  w  dw] _ 

^  Jvo  (1  I'ft 


dw] _  li^  fd^u)  ;/  d'^wj]  d^7/3 


12  Ox'^  fl2  Qg2  Q^2 


h^  r  1  O^w  d'^T/3  h^ 


12  AM  R^  Off^  Ox^  Off^  6R^ 


1  aai  ;i~2  i  aai 


O^W  0'^7]3 


0x06 0x06 


where  $  =  {(f),u,v,w)  and  'P  =  (^, »/i,?/2,?/3)-  We  remark  that  the  inner  product  for  the  state 
space  H  contains  terms  which  arise  when  considering  the  kinetic  energy  for  the  shell  while 
the  shell  contributions  in  the  V  inner  product  are  motivated  by  the  form  of  the  strain  energy 
for  a  thin  cylindrical  shell  (see  [8]).  Also,  we  note  that  with  this  choice  of  spaces  and  inner 
products,  we  can  form  the  Gelfand  triple  K  ~  with  pivot  space  H  (see  [27] 

for  a  complete  discussion  of  these  ideas). 

To  define  appropriate  sesquilinear  forms  a,  :  1/  x  V  — >  C  ,i  =  1,2,  we  group  the  stiffness 
and  wave  contributions  separately  from  damping  and  coupling  terms,  thus  leading  to  the 
definitions 


/.R  .i,\  f  ^'oh  {\0u  _  u  fOv  _  M  1  ,,  s  fdu  1  0u'\  dr/i  1  , 

^  “  -/ro  (1  -  i'^)  I  /J  ~  ^  '06  I  ^ 

,  [  Cph  \  I  \  I  Ov  w  ^  ..^“1  .1,,  ,  1  di]2\  , 

•/ro  (1  X'R  [R'^'^  R'^’'o^\  2^^  R'^\  ■^J 

f  Coh  I  1  [  1  dv  w  d«] _  h^  O^w  u  0^vi  O^rjs 

^  Jro  (1  -  ^R^  12  ^ 


h^  r  1  O^w  0^w\  0^7]3  h^ 


1  o  D2  02  i: 


12  R^  R^  06'^  dx2  06'^  6/?2 


O^W  d'^7/3 
0x06 0x06 


+  I  PS  {<i>m  —  j  dj  . 


To  account  for  the  control  contributions,  we  let  U  denote  the  Hilbert  space  containing  the 
control  inputs,  and  we  define  the  control  operator  B  €  C{U,  V*)  by 


(flu,  *)^.y  =  f  ||(/v.)„l,  ^  i  [(;v,)„|,  ^  -  1(M,)„1, 


1  r,  ^  1 

0x2  ^02  I 


10 


for  ']/  G  K,  where  y  is  the  usual  duality  pairing.  Finally,  by  letting  F  =  ^0,0.0. 

we  can  write  the  system  in  the  abstract  variational  for?n 

(yttit),  '!')  +  '!')  =  (Bu(f.)  +  F,  'i')v..,v/  ■  ( lO) 

We  reiterate  that  the  state  for  this  second-order  system  is  p(t)  =  (6(t,  •,  ■,  •),  xi{i,  •).  v{t .  ■.  ■). 
w{t,  •,  •))  in  V  H. 

In  order  to  write  the  system  in  terms  of  associated  bounded  operators,  we  first  not(' 
that  (Ti  and  (Tj  are  bounded  (there  exist  c\  and  c-i  such  that  |cri(4>,4')l  <  ci I4>|v'|'l'|\'  aiul 
|(T2(^>,^)1  <  C2|^>|v|'I'|v),  and  hence  we  can  define  operators  A\,A2  G  £(V',  V'‘)  by 

for  i=  1,2.  This  then  yields  the  abstract  second-order  system 

J/n(f)  +  A-2yi{i)  +  A\y{t)  =  Bu[t)  -f  F 


in  \/*. 

To  apply  infinite  dimensional  control  results  for  periodic  forcing  functions  to  this  problem, 
it  is  advantageous  to  write  the  system  in  first-order  form.  This  is  accomplished  by  defining 
the  product  spaces  H  =  V  x  H  and  V  =  V  x  V  with  the  norms 

and 

We  point  out  that  V  "H  ~  V*  again  forms  a  Gelfand  triple. 

The  sesquilinear  form  cr  ;  V  x  V  — >  (C  is  then  defined  by 

^(0,  X)  =  A),  (4»,  vF))  =  -  (A,  ^)y  +  .7, (T,  vF)  -f  cr2( A,  4^) 

where  x  =  'J')  and  0  =  (T,  A). 

For  the  state  y(i)  =  (y(t),yt(t))  in  ‘H-,  we  can  subsequently  write  the  system  in  the  first- 
order  variational  form 


A:)v.V  +  x)  =  {Bu[t)  -f  T{t),  x)v^y  (11) 

where  =  (0,  F(f))  and  Bu[t)  =  {0,Bu{t)).  As  usual,  the  relation  (11)  must  hold  for  all 
X  G  V.  Finally,  the  weak  form  (11)  is  /orma/Zy  equivalent  to  the  system 

yt{t)^Ayit)  +  Bu{t)  +  j^{t)  (12) 

in  Ti  where 

dom  >1=  {0  =  (T,A)  G  K  :  A  G  V,A,T  -f  A2A  G  //} 


A  = 


/ 

-A2 


(IT; 


1 


Hence  we  see  that  the  system  corresponding  to  the  original  model  can  be  written  in  various 
forms  whose  utility  depend  on  the  applications  of  interest.  For  approximating  the  dynamics  of 
the  system,  the  formulation  (7),  or  equivalently  ( 10)  or  ( 1 1 ),  is  useful  and  ultimately  leads  to  a 
first-order  matrix  system  when  considering  the  resulting  finite  dimensional  problem.  Moreover, 
as  will  be  shown  in  the  next  section,  the  weak  form  (10)  provides  a  setting  in  which  one  can 
determine  well-posedness  results  for  the  problem.  Finally,  the  first-order  infinite  dimensional 
formulation  (12)  is  an  abstract  Cauchy  equation  which  facilitates  the  application  of  infinite 
dimensional  LQR  results  to  the  problem  of  controlling  interior  acoustic  pressure  levels. 


4  Well-Posedness  of  the  System  Model 

In  this  section,  the  system  model  described  in  strong  form  in  Section  2  and  weak  form  in 
Section  3  is  shown  to  be  well-posed.  The  strategy  used  to  do  this  can  be  summarized  as  follows. 
By  employing  the  properties  of  the  sesquilinear  forms  Ui  and  (T2,  it  is  demonstrated  that  the 
operator  A  of  (13)  generates  a  Co-semigroup  S{t)  on  7i.  Due  to  the  fact  that  the  control 
and  forcing  terms  lie  in  V*  rather  than  however,  one  is  prevented  from  simply  applying 
a  variation  of  parameters  approach  to  define  mild  solutions  as  is  often  done  in  problems  of 
this  type.  Instead,  a  construction  due  to  Haraux  [17]  is  used  to  extend  the  semigroup  S{t) 
to  the  space  W*  =  [dom  A*]*  where  an  extension  of  the  variation  of  parameters  technique 
for  describing  a  mild  solution  to  the  problem  can  be  employed.  By  formulating  the  problem 
in  this  generalized  manner,  the  well-posedness  of  the  system  solution  can  be  established  in  a 
form  that  can  be  used  in  parameter  estimation  and  control  problems. 

4.1.  Generation  of  the  Co-semigroup  S{i)  on  7i 

Before  proving  that  the  operator  A  generates  a  Co-semigroup  on  7i,  we  recall  several 
definitions  that  will  be  needed  (we  note  that  these  definitions  vary  slightly  among  authors, 
e.g.  see  [24,  27],  and  that  the  definitions  given  here  will  be  used  in  the  discussion  which 
follows). 

We  say  that  the  sesquilinear  form  cr  :  K  x  1/  (L'  is  V -elliptic  if  there  exists  a  constant 

c  >  0  such  that 

Rea(<F,«I>)  >  c\^\l  for  all  <I>  €  R  . 

The  sesquilinear  form  cr  :  R  x  R  (L’  is  called  H -semiclliptic  if  there  exists  a  constant  6  >  0 
such  that 

Re(T(<D,<D)  >  6|<I>|^  for  all  <!>  e  K  . 

Finally,  we  say  that  the  sesquilinear  form  a  is  symmetric  if 

<t(<D,'I>)  =  <t('I»,<I))  for  all  vk  e  K  . 

The  following  theorem  can  be  used  to  establish  that  the  operator  A  given  in  (13)  generates 
a  Co-semigroup  on  H.  The  proof  depends  on  the  Lumer- Philips  Theorem  and  arguments  can 
be  found  in  [1]  and  pag(  s  82-84  of  [3]. 
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Theorem  1.  Let  V  ^  H  ^  H*  ^  V’  and  suppose  that  ctj  and  (T2  of  (10)  satisfy  tlu' 
properties:  <Ti  is  K-elliptic,  continuous,  symmetric,  and  (72  is  continuous,  //-semielliptic. 
Then  A  defined  in  (13)  generates  a  fVsemigroup  on  7i. 

To  see  that  our  system  (10)  satisfies  the  hypotheses  of  this  theorem,  we  first  consider  the 
sesquilinear  form  (T](^,'1')  =  (^>,'1')^  given  in  (8).  The  boundedness  of  ct]  (which  results  from 
Schwarz’s  inequality  for  inner  products  in  conjunction  with  equivalence  results  for  various 
Sobolev  norms)  was  noted  in  the  last  section.  The  V'-ellipticity  and  symmetry  of  cri  follow 
directly  from  the  definition  of  the  sesquilinear  form  as  the  V'  norm.  We  note  that  the  symmetry 
of  cr,  depends  on  the  symmetry  of  the  Donnell-Mushtari  shell  operator,  and  while  some  other 
shell  theories  such  as  that  of  Byrne,  Flugge  and  Lur’ye  also  provide  symmetric  operators 
and  resulting  sesquilinear  forms,  others  such  as  the  Love-Timoshenko  theory  do  not  yield 
symmetric  operators  (see  [19])  and  hence  would  not  fall  directly  into  our  framework  here. 
From  a  physical  perspective,  the  symmetry  of  the  shell  operator  and  resulting  sesquilinear 
form  guarantees  real  vibration  frequencies.  This  symmetry  also  guarantees  that  the  shell 
model  satisfies  the  Maxwell-Betti  Reciprocity  Theorem  which  essentially  states  that  for  a 
linearly  elastic  body,  the  general  displacement  at  a  point  m  resulting  from  a  load  at  the  point 
n  is  equal  in  magnitude  to  the  displacement  at  u  resulting  from  an  equal  magnitude  load  at 
m  (see  for  example,  [23]). 

We  next  turn  to  the  sesquilinear  form  (72  of  (9)  and  recall  that  the  boundedness  of  <72  was 
also  noted  in  the  last  section.  This  follows  from  the  definition  of  the  a2  and  the  equivalence 
of  various  norms  which  arise  from  the  components  of  the  V  norm.  The  //-semiellipticity  of 
(72  follows  directly  with  the  choice  6  =  0. 

The  results  from  Theorem  1  guarantee  that  the  operator  A  given  in  (13)  generates  a  Co- 
semigroup  S{t)  on  the  state  space  H.  Moreover,  the  semigroup  satisfies  the  exponential  bound 
^  for  t  >  0  (where  in  fact,  u;  =  0  due  to  the  fact  that  A  is  dissipative  as  shown  in 

[3]). 

In  the  case  of  a  bounded  (in  Ti)  control  input  operator,  the  usual  procedure  is  to  use  a 
variation  of  parameters  approach  to  define  a  mild  solution  to  the  system.  As  noted  in  the 
last  section,  however,  the  control  input  B  €  C{IJ,V')  defines  the  product  space  control  term 
Bu(t)  =  {0,Bu{t))  G  {0}  X  V*  <Z  V  V*  =  V*.  Since  Bu{t)  lies  in  V’  rather  than  in  H.  the 
usual  variation  of  parameters  ideas  are  not  feasible.  We  are  therefore  motivated  to  extend  the 
semigroup  S{t)  on  to  a  semigroup  S{t)  on  a  larger  space  W*  D  {0}  X  V'*  so  that  the  mild 
solution 

is  well-defined  for  Bu-j-F  G  L^((0,  T),  K*)  (in  the  following  development,  it  will  be  shown  that 
the  space  W*  is  given  by  W*  =  [doni-4*]*).  This  then  provides  a  setting  in  which  to  guarantee 
the  well-posedness  (including  continuous  dependence  on  initial  data  and  non  homogeneous 
input  terms)  of  the  solution. 


4.2.  Description  of  and  dom  A‘ 


Before  discussing  the  extension  of  S{t)  from  7i  to  W*,  it  is  useful  to  describe  A“  and 
dom  A*.  We  first  recall  from  (13)  that  A  and  dom  A  are  given  by 

dom>l={0  =  (T,A)eW:  Ae  +  H} 

(  -^,T-A,A  )  • 

From  the  usual  definition  of  the  adjoint  we  now  want  to  find  dom  A*  and  A’'  satisfying 
dom  >1*  =  {x  =  (^,  €  H\^  e-v  (^0)(x)  is  continuous  on  dom  A} 

and 

Me,x)„  =  (eM-x)« 

for  all  0  =  (T,  A)  e  dom  A  and  x  =  '!')  €  dom  A*.  It  follows  directly  from  the  definition 

of  the  adjoint  operator  that  x  G  dom  A”  if  and  only  if  there  exists  T  =  (71,72)  ^  H  =  V  x  H 
such  that 

(^e,x)H  =  (0.r)K 

for  all  0  €  dom  A.  The  expansion  of  this  relation  yields  the  condition 

(A,<I>)v.  +  {{-A,r  -  A2A),^)m  =  (T,7i>v  +  (A,72)h  (1-^^) 

for  all  0  =  (T,  A)  €  dom  A  and  x  =  (^,'^)  €  dom  A*.  By  noting  that  (15)  must  hold  for 
A  =  0,  it  follows  that 

{-/1,T,4)„  =  (T,7,)v  =  {/1,T,7,)„ 

for  all  T  e  dom  Ax.  The  second  inequality  results  from  the  observation  that  (T,7i)y  = 
(Ti(T,7i)  =  (/llT,7l)^/.  y  =  (AiT,7i)^  for  0  =  (T,0)  e  dom .4.  It  follows  immediately  that 
(AiT,  — ^  ~  1\)h  ^  ^1’  which  implies  that 

7,  =  -'F 

due  to  the  K-ellipticity  of  <Ti.  Moreover,  71  G  V  implies  that  'F  G  V. 

When  7i  =  is  substituted  into  (15),  the  relationship 

(/liA,<I>)^  +  (-A2A,'I')„  =  (A,72)„ 

results  for  all  A  G  V'.  From  the  last  equality,  it  follows  that  the  mapping  A  (4iA,4>)^  + 
{—A2A,'^)ff  is  continous  in  H.  This  in  turn  implies  that  —  4^^  G  H  and  that 

(A,/i;4> -/!;♦>„  =  (A, 7j)„ 

for  all  A  G  K.  Hence  we  have 

72  =  4t<I>  -  Al<if  . 
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Thus  A'  and  doni  A*  are  given  by 


dom  A''  =  {x  =  (<I>, '!')  e  K(4'  G  V,  -  /I*'!'  e  //} 


-4*X  = 


-vp 


(16) 


-  A'^^  } 

Finally,  we  note  that  because  'I'  6  F  for  any  x  ^  dom  A’,  it  follows  immediately  that 

dom  A*  C  V  X  V  =  V  . 


4.3.  Extension  of  S{t)  to  [dom  .A*]* 


With  A’  and  dom  A’  described  as  above,  we  now  <'onsider  the  extension  of  the  semigroup 
S{t)  generated  by  A  from  Ff  to  a  larger  space  W*  which  is  defined  below.  This  will  be 
accomplished  by  employing  the  extrapolation  space  techniques  discussed  by  Haraux  in  [17]. 
Other,  more  detailed  examples  using  these  techniques  can  be  found  in  [5,  6]. 

The  space  W  =  [dom  .4’]  is  taken  to  be  dom  A*  with  the  inner  product 


for  some  arbitrary  but  fixed  Aq  with  Xq  >  uj  (recall  that  the  original  solution  semigroup  satisfies 
the  bound  j<5(<)|  <  As  proven  in  [6],  the  resulting  W  norm  is  equivalent  to  the  graph 

norm  corresponding  to  that  is,  there  exist  constants  c\  and  such  that  for  any  4>  €  W, 
we  have 


(If) 


l<I>lw  +  |>l*<f|K<c.2|(Ao-.4*)<I>|„  . 


As  discussed  in  [5,  6],  the  space  W  defined  in  this  manner  and  with  this  norm  is  densely  and 
continuously  embedded  in  Ti.  Hence  we  can  formulate  the  Gelfand  triple  W  ■— >  ~  K*  >  W* 
with  pivot  space  7{.  Moreover,  the  dual  space  W*  is  isomorphic  to  the  completion  of  with 
respect  to  the  norm  |A|w*  =  |(Ao  —  A)~^h\ji. 

To  facilitate  the  arguments  for  extending  the  operator  A  from  dom  A  QH.  io  all  of  ,  we 
define  the  sesquilinear  form  cr  :  Ff  x  W  €  by 


a(0,x)s(0,>l*x)„ 

for  all  0  G  X  ^  Due  to  the  equivalence  of  the  W  norm  and  the  graph  norm  corre¬ 
sponding  to  A*.,  it  follows  that 

|d’(0,x)l  <  |€)|kM’x|w  <  <’|G|w|x|i^  • 

This  then  implies  that  for  each  0  €  Ti,  the  mapping  x  ^  d’(0,  \)  is  in  W”.  Hence  we  can 
define  Ai0  G  W’  by  (.40^  (x)  =  d‘(0,x)  or  equivalently, 

(.4e)(o  =  (e,>C\)„ 


15 


for  all  0  G  X  G  W.  From  this  definition,  it  follows  that  A  G  C{'H,W*)  with  dom  A  =  1-1. 
Moreover,  if  0  G  dom  it  follows  from  the  Riesz  representation  theorem  that 

(ie)  (v)  =  =  {Ae.x)n  =  (^0)(\) 

for  all  \  G  dom  A".  Thus  we  see  that  A  is  an  extension  of  the  original  operator  A  from 
dom  A  <Z  'H  io  aW  o{  l-L. 

Having  extended  the  operator  to  the  full  state  spare,  our  final  objective  in  this  subsection 
is  to  observe  that  A  is  the  infinitesimal  generator  of  a  f'o-semigroup  on  W*.  The  proof  of 
this  result  is  given  in  [6],  and  we  simply  state  the  conclusion  in  the  following  theorem  while 
directing  the  reader  to  that  paper  for  details. 

Theorem  2.  The  operator  A  is  the  infinitesimal  generator  of  a  T’o-semigroup  S{t)  on  W” 
which  is  an  extension  of  S{t)  from  7^  to  W*. 


Having  described  how  the  original  semigroup  S{i)  can  be  extended  to  a  semigroup  S{t)  on 
VV*  =  [doniwA*]*,  we  now  want  to  show  that  the  space  {0}  x  V’  which  contains  the  control  and 
force  input  terms  is  itself  contained  in  VV*.  It  is  tempting  at  this  point  to  recall  that  we  have 
the  Gelfand  triples  V  Tf  ~  'H*  ^  V*  and  VV  ~  Ti*  »  VV’  as  well  as  the  embedding 
VV  C  V,  and  from  this  assume  that  the  desired  result  is  obvious.  We  point  out,  however,  that 
while  VV  C  V,  the  embedding  is  not  in  general  continuous  which  implies  that  VV  >  V  can  not 
be  assumed  to  be  true  (if  A*  is  V-elliptic  for  example,  one  can  obtain  the  desired  continuity 
of  the  embedding;  however,  this  is  not  the  case  in  this  problem).  Hence  more  care  must  be 
taken  in  order  to  argue  that  W*  D  {0}  x  V* 

To  this  end,  let  h  =  (0,  A*)  be  an  arbitrary  functional  in  {0}  x  V*.  In  order  to  show  that 
h  G  VV*,  we  need  to  show  that  y  =  (^>,  'V)  /t(x)  is  bounded  on  VV.  Now,  it  follows  that 

\h{x)\  =  |A*('I')1  <  lA’Iv-lvVlv  <  mv  .  (18) 

Moreover,  from  the  definition  (16)  of  A*,  we  see  that  A'x  =  (  — which  implies 
that 

I'I'lv  <  M'xlw  S  +  I.X'Ih  <  C2|x|w  •  (19) 

(The  latter  inequality  results  from  the  equivalence  of  the  VV  norm  and  the  graph  norm  cor¬ 
responding  to  A*  as  described  in  (17).)  From  (18)  and  (19)  it  follows  that  there  exists  a 
constant  c  such  that 

IMx)l  <  c|x|w 

therefore  showing  that  h  G  W*.  We  have  thus  established  the  following  result  which  is  essential 
in  developing  expression  (14). 

Lemma  1.  Under  the  assumption  of  Theorem  1,  we  have  {0}  x  V'*  C  W*  =  [dom  >4*]*. 

In  summary,  we  are  now  able  to  show  the  following.  The  system  (11)  is  well-posed  in  the 
sense  that  (14)  is  a  well-defined  entity  that  can  be  taken  as  the  mild  solution  to  (11).  More¬ 
over,  because  (14)  is  well-defined,  the  usual  theorems  for  continuous  dependence  with  respect 
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to  initial  data  and  nonhoniogeneous  terms  (control  and  external  forces)  follow  ininiediately. 
This  facilitates  the  discussion  of  approximation  and  control  ideas  for  these  problems  (sc'e  for 
example,  [3,  4]). 

In  the  event  that  cr-i  of  (10)  is  also  V'-elliptic,  one  can  argue  (see  [5]  for  details)  existence  of 
a  weak  or  variational  solution  to  the  system  (10)  which  agrees  with  the  mild  solution  obtaiiual 
from  (14).  In  the  exa?nple  of  Section  3,  the  form  <72  is  V'-elliptic;  nevertheless,  one  can 
turn  to  standard  results  from  the  theory  of  Co-semigroups,  [22.  26],  to  obtain  existenc<'  and 
uniqueness  of  solutions  for  the  system  under  consideration  in  this  paper.  For  example,  we 
have  the  following  result. 

Theorem  3.  Ckinsider  the  system  represented  by  (10),  (12)  or  ( 14)  and  suppose  that  the  map¬ 
pings  t  H- >  u[t)  and  t  1-^  F{t)  from  [0, T]  to  IR’  and  V’',  respectively,  are  Lipschitz 
continuotis.  Then  for  each  To  ^  H  =  dom  A,  we  have  that  (12)  taken  with  T'lO)  =  To 
has  a  unique  strong  solution  given  by  (14). 

This  theorem  is  readily  established  by  appealing  to  Corollary  2.11,  p.  109  of  [22].  in  the* 
context  of  our  arguments  above.  Letting  G(t)  ~  (0,  Bu{t)  -f  F{t))^  and  considering  S{t)  on 
the  reflexive  Banach  space  .V  =  W”,  we  can  argue  that  (i  is  Lipschitz  continuous  under  the 
hypotheses  of  the  theorem.  In  particular,  the  inequalities  (18)  and  (19)  inqily 

|C(<)|w  <  c2\Bu{t)  +  F{t)\v  , 

so  that 

|6’(0  -  C(.s)|w  <  C2  {lei  W)  l«(0  -  «(•'*)!  +  1^(0  -  • 

It  follows  that  (14)  provides  the  strong  solution  (i.e.,  differentiable  a.e.  in  the  W*  sense  - 
see  [22])  to  (12)  interpreted  in  the  W'  sense. 


5  Conclusion 

In  this  paper,  a  modeling  set  of  partial  differential  equations  de-scrilung  tin*  dynamic.s  (d  a 
structural  acoustics  problem  was  presented,  and  a  mathematical  framework  amenable  to  the 
development  of  approximation  schemes  for  forward  simulations.  ])arameter  ('stimation.  and  the 
application  of  PDE-based  feedback  control  strategies  was  tlevelop«'d.  The  stru<  ture  consists 
of  a  thin  cylindrical  shell  with  hard  caps  at  the  ends,  and  this  component  of  the  system  was 
modeled  by  the  Donnell-Mushtari  shell  equations.  A  constitiitive  law  postulating  that  strt'ss 
is  proportional  to  a  linear  combination  of  strain  and  strain  rate  was  assum<*d  which  yit'lds  a 
Kelvin- Voigt  type  of  damping  in  the  shell.  The  structural  dynami<  s  were  tln-n  coupled  to  the 
interior  acoustic  field  through  pressure  and  momentum  balance  conditions.  Finally,  this  nnxlel 
includes  the  contributions  due  to  the  activation  of  piezoceramic  patches  which  are  I)unded  to 
the  shell  and  ultimately  will  be  used  to  control  the  interi(»r  acoustic  pn'ssun*  h-vels. 

To  accommodate  the  presence  of  the  |)iezocerami<'  patches  and  difft-ring  mati'i  ial  i)roi)ei  i  ies 
of  these  patches,  the  material  parameters  for  the  combined  structun'  should  lx*  taken  to  be 
piecewise  constants.  This  leads,  however,  to  difficulties  in  the  strong  form  t)f  the  system 
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equations  since  these  parameters  are  contained  in  moment  and  force  resultants  which  are 
differentiated  when  forming  the  equations  of  motion  for  the  shell.  Moreover,  the  discontinuities 
introduced  by  the  patches  lead  to  an  unbounded  control  operator  since  it  involves  derivatives 
of  the  Heaviside  function  and  Dirac  delta.  To  avoid  these  problems,  the  weak  form  of  the 
system  equations  was  also  presented. 

In  the  weak  form,  the  derivatives  appear  on  the  test  functions  instead  of  on  the  moments, 
thus  eliminating  the  difficulties  associated  with  the  discontinuous  parameters  and  patch  inputs. 
The  weak  form  is  also  advantageous  for  many  approximation  schemes  since  it  reduces  the 
smoothness  requirements  for  the  basis  elements.  Finally,  when  the  weak  form  is  posed  in  the 
context  of  sesquilinear  forms,  convergence  and  well-posedness  issues  can  be  considered. 

The  first  step  in  proving  the  well-posedness  of  the  system  model  was  to  verify  that  the 
first  order  system  operator  A  :  dom  A  C  'H  Ti  generatt'd  a  fVsemigroup  S{t)  on  Ti.  This 
semigroup  was  then  extended  via  the  extrapolation  techniques  described  by  Haraux  to  a  larger 
space  W*  =  [dom  A']"'  so  as  to  be  compatible  with  the  force  and  control  inputs  which  lie  in 
{0}  X  V’.  Finally,  the  mild  solution  in  terms  of  the  extended  semigroup  was  formulated  and 
existence  and  uniqueness  results  were  obtained.  Hence  the  framework  and  model  presented 
here  is  mathematically  well-posed  cis  well  as  amenable  to  the  development  of  computational 
strategies. 
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